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0\ . Abstract 



The impact ionization rate and its orientation dependence in k space is calculated for ZnS. The 



s 

' numerical results indicate a strong correlation to the band structure. The use of a q-dependent screening 

O , function for the Coulomb interaction between conduction and valence electrons is found to be essential. 



A simple fit formula is presented for easy calculation of the energy dependent transition rate. 

PACS-No. 72.20.Ht, 72.20.Jv 

The impact ionization rate is of basic interest in understanding the electron transport in semiconductors 
at high electric field strengths. This microscopic quantity is necessary for calculating the probability of the 
electron to ionize which is described by the ionization coefficient a. Interesting features are the orientation 
dependence and the threshold behavior, which were studied for Si |0, ||, ^, Q and GaAs ||. 

Recently, we have determined the ionization rate for GaAs and Si following a scheme presented by 
Sano and Yoshii Q for the evaluation of the impact ionization rate for Si using a very efficient deterministic 
numerical procedure. Their special integration scheme takes into account a rather high number of points in 
the Brillouin zone so that they were able to study effects of the anisotropy in k space as well as the threshold 



behavior for Si. Both were found to be direct consequences of the band structure. 

In this paper we investigate the impact ionization rate and its orientation dependence in the wide band 
gap semiconductor ZnS {Eg = 3.63 eV). The calculation is performed similar to [|| utilizing the integration 
method of Sano and Yoshii again. 

Inserting a band structure which is obtained from local pseudopotential parameters for ZnS, the cor- 
relation of the particular band shapes with the impact ionization rate is shown. Exchange and umklapp 
processes are included in the evaluation of the impact ionization rate. The consideration of an appropriate 
wave-vector dependent screening function for the Coulomb interaction in semiconductors is of significant 
influence on the results. 

The impact ionization rate is usually determined from Fermi's golden rule, 

x5[E,,{\^^) + E,,{\^2) - E,,{\^^) ~ E,M] , (1) 

where £7 and Vi arc the crystal volume and the band index, resp. Mtot denotes the matrix element for 
ionization including direct (D), exchange (E) and umklapp processes, 

iMtotP = 2|Afo|' + 2|Mb|2 - |Af* Mb + Af^M* | . (2) 

The initial electron states belong to the conduction (1) and valence band (2), whereas the final ones are 
both in the conduction band [(3) and (4)]. The integrals are extended over the Brillouin zone. The ionization 
rate, Eq. (|^), is integrated over all directions in k space to get an averaged, only energy dependent impact 
ionization rate 

E/d^ki 5[E,,{\^i)-E] r(ki,iyi) 

^^^^ " E/rf^ki <5[i?.,(ki)-£;] ■ 

The analytical evaluation of the transition rate |7| requires strongly simplifying assumptions such as 
a parabolic band structure and free electron wave functions, which are not well founded for semiconductors 
with a direct (GaAs, ZnS) or even an indirect gap (Si, ZnS when omitting the first conduction band which 
is not contributing). Therefore, we numerically evaluate Eqs. (0) and (^) in order to take into account a 
realistic band structure and corresponding wave functions in the entire Brillouin zone. This is done using 



a grid method following Sano and Yoshii (Ref. 2), who made extensive use of symmetry relations imposed 
by the crystal structure and were able to take into account a great number of points in the Brillouin zone 
for the evaluation of the integrations in Eq. (|l|) for Si. This method is accurate enough to treat even more 
sensitive effects such as the threshold behavior or the orientation dependence of the interband transition rate 
successfully. 

We have calculated the band structure and wave functions within the frame of the empirical pseudopo- 
tential method H], utilizing the set of pseudopotential parameters from Ref. 9. The expansion runs over 
113 G vectors to produce a reasonable band structure in the entire Brillouin zone and, especially, to give a 
realistic description near the threshold. A plot of the band structure is presented in Fig. 1. 

The interaction between the electrons in the bands is described by a statically screened Coulomb potential. 
The consideration of the wave-vector dependence of the dielectric function e{q) is of importance for the 
calculation of the matrix elements in semiconductors as pointed out already by Laks et al. For the 

interaction between conduction and valence electrons, we apply a q-dependent dielectric function derived 
by Levine and Louie | pT[ . This model incorporates the correct long-range and short-range properties of the 
dielectric function and, in addition, some consequences of the presence of a gap in the excitation spectrum as 
characteristic of semiconductors. The interaction between electrons in the conduction band can be modeled 
by a Debye potential with an inverse screening length n — [noe^ / {eokBT)]^^^^ for a temperature of T = 300 
K and a corresponding electron density of tiq = 10^^ cm^'^. The choice of uq is of less importance, as there 
are very few ionization events with small momentum transfers oi q ~ k. 

The k vectors are restricted to the first Brillouin zone with respect to the reduced zone scheme. We 
want to mention at this point that for umklapp processes we have to consider terms which go beyond 
the expansion up to the 113 G vectors of the pseudopotential band structure calculation. The respective 
expansion coefficients, i.e. those belonging to reciprocal lattice vectors with |G| > V20{^), were calculated 
by means of perturbation theory |l2j . 

For the numerical treatment of Eqs. (|l|) and the (5-function is replaced here by rectangles of unit area 
with a height of 1/SE and a corresponding width of SE. A value of 6E = 0.1 eV was found to be large 
enough to ensure good convergence. This was proved by test calculations with a finer grid as well as with a 
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value of SE — 0.2 eV. Of course, the broadening factor is related to effects such as collision broadening and 
intra-coUisional field effect However, the study of such effects is not intended here. 

An important feature of the grid method used in our previous paper ||] and in Ref. 2 is, that the storage of 
energy values and wave functions can be restricted to the irreducible wedge (IW) if the symmetry properties 
of the Brillouin zone are fully exploited. The wave function for a wave vector Tkk, which is obtained by 
operating one of the 48 transformations of the lattice point group on a wave vector k of the IW, is then 
determined from ^^(k) (see [§, D). The mesh in k space is constructed with regard to the limits of the IW, 

277 

< k^ + ky + k,< l-5(— ) , 

< k,<ky<k,<{ — ), 

as well as to the condition for the interval length of the assumed grid Afc^ = ^(^), i — {x, y, z} with n e H. 

For the determination of the ionization rate r(k, J/), we take a mesh of 152 points in the I W [Afc^ = 
For the study of the average ionization rate R{E), the summations over the conduction band indices were 
restricted to the four lowest bands, whereas all four valence bands were taken into account. 

In order to investigate the influence of the band structure (see Fig. 1) on the ionization rate, we examine 
r(ki,j/i) [Eq. (|l|)], especially along three symmetry lines. In the first conduction band there are no points 
at all satisfying the threshold condition and, therefore, the total rate is due to electrons in the higher - 
mainly the third and fourth - conduction bands. The anisotropy of the ionization rate of ZnS for electrons 
initiating from the second, third and fourth band is shown in Fig. 2. Comparing these results with the 
band structure data, a coincidence of the qualitative behavior of the impact ionization rate for the three 
directions and the third and fourth conduction bands can be stated. However, the attempt to manifest this 
correlation by applying a power-law fit formula for r(ki, vi) as used in our prevoius paper |^ did not lead to 
a satisfactory result. The validity of this Keldysh-type fit strongly relyes on the isotropy of the matrix 
elements. Contrary to calculations for GaAs ^ and Si |^ ^, these exhibit a notable dependence on the k 
value of the initial electron, which hence does not allow for such a fit. 

Furthermore, the necessity to consider the region of enlarged energies above the second conduction band 
where band degeneracies and multiple level crossings occur, complicates the correlation of impact ionization 
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to the band structure. For such high energies, also the quahty of empirical pseudopotential band structures 
becomes questionable. The desirable application of a nonlocal band structure may lead to more rigorous 
results for the impact ionization rate. 

The average energy dependent ionization rate R{E) can, however, be approximated by such a simple fit 
formula, which is given by 

R{E) = P[E - EthT , (4) 

where the parameters are fixed to P = 5.14 x IQ^^^ s^^{eV)~°' and a = 5.183. The threshold energy was 
set to Eth — 3.8 eV (The calculation gave first ionization events above Eth — 3.93 eV; regarding Eth as a 
fit-parameter, best fits may be obtained with Eth ~ 3.5 eV). Fig. 3 shows the result as obtained from our 
numerical calculations together with the fit R{E). The increase of the curve indicates a threshold behavior 
similar to that obtained for GaAs, a typical wide band-gap semiconductor with a hard threshold. 

We have neglected in our present calculations phonon assisted impact ionization processes as well as deep 
level ionization which may become of importance at high fields. Furthermore, other high-field effects such 
as collision broadening or the intra-collisional field effect have to be considered for a more comprehensive 
description of the transport process of hot electrons in semiconductors similar to [|l3[ |l5) . 

The aim of further calculations may be the study of the influence of the obtained anisotropy effects on 
the ionization coefficient, where an additional effect of the band structure is expected Furthermore, the 
influence of dynamic screening which was found to be substantial for Si [0 has to be studied also for ZnS. 
The inclusion of the obtained impact ionization rate would also improve present calculations for high-field 



transport in ZnS |18, hg 
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FIGURE CAPTIONS 

FIG. 1: Pseudopotential band structure of ZnS. The labels indicate the conduction bands including degen- 
eracies. 

FIG. 2: Ionization rate for ZnS calculated for k points along three symmetry lines. 

FIG. 3: Average ionization rate for ZnS obtained from the present numerical calculation and from the fit 
formula (Eq. indicating the contributions of the lowest four conduction bands. 
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